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We study dynamics of filling of an initially empty finite medium by diffusing particles A and B, which 
arise on the surface upon dissociation of AB molecules, impinging on it with a fixed flux density /, 
and desorb from it by the reaction A + B — > AB — > 0. We show that once the bulk diffusivities differ 
(p = Da/ Db < 1), there exists a critical flux density I c (p), above which the relaxation dynamics 
to the steady state is qualitatively changed: on time dependences of CAs/c e (c e being the steady 
state concentration at t — > oo) a maximum appears, the amplitude of which grows both with I and 
with Db/Da ratio. In the diffusion-controlled limit / > / c at p < 1 the reaction "selects" the 
universal laws for the particles number growth Ma = Mb oc i 1 / 4 and the evolution of the surface 
concentrations ca s oc i -1 '' 4 , cb s oc i 1//4 , which are approached by one of the two characteristic 
regimes i) and ii) with the corresponding hierarhy of the intermediate power-law asymptotics. In 
the first of these cas goes through a comparatively sharp max(cA s /c e ) oc 7 1//6 , the amplitude of 
which is p-independent, in the second one cas goes through a plateau-like max(cAs /c e ) ocp- 1/4 , the 
amplitude of which is /-independent. We demonstrate that on the main filling stage the evolution 
of the M(t)/M e ,CAs(t) /c e , and cg s (t)/c e trajectories with changing p or J between the limiting 
regimes i) and ii) is unambigouosly defined by the value of the scaling parameter fC — p :i ^ 2 J (J 
being the reduced flux density) and is described by the set of scaling laws, which we study in detail 
analytically and numerically. In conclusion, we analyze specific features of the long-time relaxation 
dynamics and calculate exactly the relaxation rate u>(p,J). 
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I. INTRODUCTION 

For the last decade the reaction-diffusion system A + 
B — > 0, where unlike species A and B diffuse and irre- 
versibly react in the bulk of a d-dimensional substrate, 
has acquired the status of one of the most popular objects 
in nonequilibrium statistical physics (!]]. Two situations 
have been investigated most intensively: (i) dynamics of 
the A + B — > annihilation in an infinite system with 
initially homogeneously (randomly) and equimolarly dis- 
tributed reactants, in which below the critical dimen- 
sion d c = 4 a dynamical clustering develops (fluctuation- 
induced like-particle domains formation) and, as a result, 
an anomalous reaction deceleration arises (Ovchinnikov 
and Zeldovich, 1978 @, Toussaint and Wilczek, 1983 @); 
(ii) behaviour of the dynamic reaction front in an infinite 
system with initially spatially separated reactants (Galfi 
and Racz, 1988 Q), and structure of the steady state 
front in a finite system, at the ends of which are injected 
equal currents of A's and B's particles (Ben-Nairn and 
Redner, 1992 Cornell and Droz, 1993 §). 

Recently, we have shown j^U, that m another wide 
class of RD system, where reaction and diffusion are spa- 
tially separated ( i.e. reaction proceeds on the surface 
of the medium and diffusion proceeds in its bulk ) the 
interplay between reaction and diffusion acquires qual- 
itatively new features and leads to the threshold self- 
organizing dynamics of the A + B — > 0. It has been found 
that once particles A and B diffuse at different mobilities 
from the bulk of finite medium onto the surface and die 
on it by the reaction A + B — > 0, there should exist some 



threshold difference in the initial numbers of A and B 
particles , A c , above which the process of their death , 
instead of usual deceleration, starts to accelerate autocat- 
alytically. Moreover, it has been demonstrated that in 
the diffusion-controlled limit A — * oo a new critical phe- 
nomenon develops in the system - annihilation catastro- 
phe, which arises as a result of self-organizing explosive 
growth (drop) of the surface concentrations of, respec- 
tively, slow and fast particles (concentration explosion) 
and manifests itself in the form of an abrupt singular 
jump in the desorption flux relaxation rate. In the limit 
of strong difference of diffusivities the annihilation catas- 
trophe leads to the phenomenon of abrupt disappearance 
of the flux (flux breaking effect) , which may pretend to be 
one of the most dramatic manifestations of the reaction- 
diffusion interplay. 

The aim of this paper is, in contrast to the pure annihi- 
lation problem A + B — > 0, to consider dynamics of filling 
of an initially empty finite medium by diffusing particles 
A and B, which arise on the surface upon dissociation of 
AB molecules, impinging on it with a fixed flux density /, 
and desorb from it by the reaction A+B — > AB — > 0. We 
assume that the density of the input flux / is not too large 
and (or) the reaction rate constant is sufficiently large so 
that in the process of filling up to the steady state the sur- 
face coverages by adatoms A a d s and B ac [ s remain small 
enough and, therefore, the reflection of AB molecules 
from the occupied sites can be neglected. In addition, 
we assume that the particles exchange rates between the 
surface and the subsurface layer are great, and already at 
early stages of diffusion into the bulk a quasiequilibrium 
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is reached between the surface and the subsurface layer. 
Finally, we assume that in our 2d-reaction + 3d-diffusion 
problem diffusion smoothes out lateral fluctuations and, 
hence, planar distribution of A's and B's particles is sus- 
tained mesoscopically uniform, so that 1) the desorption 
flux density may be described by the mean-field expres- 
sion and 2) the problem can effectively be considered as 
one dimensional. In the framework of these assumptions 
we formulate a closed system of nonlinear boundary value 
diffusion problems, which we then investigate in detail 
analytically and numerically. 

II. MODEL 

Let the both surfaces X = ±£ of an initially empty in- 
finitely extended slab of thickness 2£ be hit with a fixed- 
density I flux of AB molecules, which dissociate on unoc- 
cupied sites to adatoms A a ds and Bads with probability 
f. Adatoms A a d s and B a d s migrate along the surface 
and either, on running into one another, irreversibly des- 
orb by the reaction A a d s + B a d s — > AB — > or diffuse 
into the bulk, gradually filling it to some steady state 
concentration c e (I) ( Fig.l ). Because of planar spatial 
homogeneity (see below), the bulk diffusion flux must 
be directed normally to the surface plane, that is, the 
problem is effectiveliy one dimensional. The boundary 
conditions for the bulk diffusion equations 

d t Ci = Did^ct (1) 

(i = A, B) can be derived from the balance of the flux 
densities on the surface 

Pi = I a -I d - If (2) 

and in the subsurface layer 

ac ts =I? b -lP\ s . (3) 

Here pi is the surface concentration of i-adatoms 
(cm -2 ), Ci S = Ci \x=e is the concentration of i-particles 
in the subsurface layer (cto~ 3 ), a is the lattice parameter, 
I a and I d are adsorption and desorption flux densities, 
I? b — I?—I b — rfpj— r|cj s ^o is the surface-to-subsurface 
layer flux density, where r? and T b are the rate constants 
of i- particles transition from the surface into the subsur- 
face layer and back, respectively, 9 = 1 — 9 a — Ob is the 
fraction of vacant sites on the surface, which we take to be 
close to unity {Oi <C 1), if \ s — DidxCi \x=i is the diffu- 
sion flux density at the surface (by symmetry, we consider 
the interval [0,£] only with the condition dxCi |x=o= 0). 
We assume that the subsurface layer - surface barrier is 
not much different from the diffusion barrier in the bulk 
(T b ~ Di/a), and probability of the adatoms transition 
to the subsurface layer is much greater than probabil- 
ity of their desorption (Z^ ^> I d ). Then, prior to des- 
orbing, the particles multiply go to the subsurface layer 
and back, and already at early stages of diffusion into 



the bulk a quasiequilibrium must be established between 
the surface and the subsurface layer, pi ~ /jCj s , where 
fi = T b /Tf is the surface segregation coefficient. More- 
over, if the surface segregation is not enough strong (at 
elevated temperatures usually / ~ (1 — f0 2 )a), and the 
system size is quite large, £ 3> fi, then at comparatively 
short (in £ 2 /Di scale) times fi < if = < t, 

when the number of particles in the bulk (per surface 

unit) Mi = J Q CidX ~ Ci s \J D{t much exceeds that on 
the surface /iC, s , the capacities of the subsurface layer 
(if \„= If) and of the surface {fi\ci s \ < If \ s ) may be 
neglected, therefore conditions (2) and (3) are degener- 
ated to the following 

Mi - if \.= r - i d . 

According to Refs. [fTo|| , Q| , and 0], in a distributed 
d- dimensional system A + B — > with injection of A and 
B particles, with their difference being strictly conserved 
(correlated landing), a critical (marginal) dimensional- 
ity, above which no fluctuation-induced segregation of A 
and B occurs, equals d c = 2 (2c?-reaction + 2<i-diffusion). 
ft is to be expected that in our case (2d-reaction + 3d- 
diffusion) 3d-diffusion smoothes out lateral fluctuations, 
and a planar distribution of A and B particles is sus- 
tained mesoscopically uniform. So, the desorption flux 
density may be described by the mean-field expression 
I d = kpAPB- From the condition If b = I a — I d it follows 

Pl = hc ls e [\ + {i a -i d )/i% 

therefore, taking I a — I d <C I b , we obtain 

i d = Kc AsCBs el[\ + (r - i d )(i/i b A + 1/4)], 

where k = kf^fs is the effective reaction rate constant 
(cm 4 /s). The above requirement I a ,I d <C X? ~ X b im- 
poses on the concentrations limitations I a /T b A <C ca s *C 
r^/K,/ a /r^ < c Bs < T b A /n, the lower boundary of 
which defines the conditions for the crossover to the 
subsurface layer-surface quasiequilibrium regime. Re- 
strictions on the flux densities I d ,I a <C T b A T b B /K follow 
herefrom. Taking into account the reflection from oc- 
cupied sites, we write I a — IB\ to obtain I a — I d = 
x0q{I - kc As c Bs ), where factor x = [1 + 9 ((n/T b B )c A s + 
(ft/r^GBs)] -1 ~ 1. So, after a short-term transient 
stage t ^S> t tr — max(/ 2 /Di) the boundary conditions 
take the form 

Mi = Did x ci \ s = x0l(I - kcasCbs), (4) 

where 1 — x@o "C 1. We assume that on a transient 
stage the desorption can be neglected, i.e. t tr <C t q , 
where t q is a characteristic time at which the desorption 
flux becomes comparable with the input flux I d ~ /. 
Using then the Laplace transform, at t <C t q , tf — £ 2 /Di 
one can easily obtain the complete solution of (l)-(3), 
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wherefrom at t » Uf >• t is (here t is = D i /(T b i ) 2 and 
t if = ff/Di) it follows 



M=M (0) 
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where Afj® — It and c 



(o) 



H= are the solutions of 

\J-KUi 

(1),(4) with the initial conditions Ci(X, 0) = at x#o = 1 
and I d /I — > 0. From (5) it is seen that the influence 
of the transient stage rapidly decays in time. So, as 
we are mainly interested here in the system's behavior 
at f > tt r , we shall take Ci(X,0) = as the initial 
condition for (1), (4). According to (5), the condition 
I d <~ nl 2 t/ 1 \JTJaDb <C / is reduced to the requirement 

t<-t q = VD i a k Db , whence it follows I < ^/D A D B /Kt tr . 
Introducing the index "7?" (heavy) for the slower diffus- 
ing species and the index "L" (light) for the faster one, 
taking x#o = 1 an d going to dimensionless variables, we 
come finally to the boundary value diffusion problem 



dh/dr = V 2 h , dl/dr= (l/p)V 2 Z, 



V/i |.= (l/p)V*| a =J-W. 



(6a) 



(6b) 



with the conditions of symmetry V(/i, I) \ x =o= 
and the initial conditions h(x, 0) = Z(x, 0) = 0. Here 
h(x,r) — Ch/c* and 1(x,t) = Cl/c* are the reduced 
concentrations, J = J/J* is the reduced flux density, 
V == d/dx x = X/£ e [0,1] is the nondimensional 
coordinate, r = Dnt/£ 2 is the nondimensional time, 
p = Dh/Dl < 1 is the ratio of species diffusivities, 
h = kc 2 = (Dn/Vkl) 2 and c* = Dh/k£ are the charac- 
teristic flux density and concentration scales, from which 
the relaxation dynamics to the steady-state becomes 
diffusion-controlled (in the reaction-controlled regime the 
both species are distributed uniformly ci s = M /£ = It/ 1, 
so, from the condition I d <~ n(It/t} 2 <C I it follows 
tn ~ £/VkJ. By comparing tn with the characteristic 
diffusion time of heavy species tf ~ £ 2 /Dh, we obtain 
J* = (Dnl^l) 2 ). 

According to (5), with an accuracy to a neglibly small 
capacity of the surface layer (r 3> r tr = max (d 2 H , pd 2 L ) , 
di = fi/£ <C 1), equal amounts (per surface unit) of H 
and L particles diffuse into the bulk Mr = Al = A/", 
therefore in the steady state (r — > oo) 



h, 



l e = N e = \fj, 



where N = Af '/AT* is the reduced number of parti- 
cles and AC = c*l = Dh/k is the characteristic scale 
of the number of particles. In dimensionless variables 
requirement 1 — \ "C 1 with account of Y\ ~ Di/a 
leads to the conditions h s <C T b L £/Dn <~ p _1 (^/a) and 
Z s <§C T b H £/DH ~ £/a whence there follows a limitation 
on the reduced flux density Z e = h e = \fl <C £/a and, 
therefore, J € ~ (£/ a ) 2 - The requirement 0, <~ 
a 2 pi <C 1 leads to the conditions h s <C (n/a 2 DH)d H 1 
and Z s <C (k / 'a 2 'D 'jj)rf^ , whence by assuming <iff ~ dz, 
we obtain J « J 8 " ~ (k / a 2 D H ) 2 d H 2 . Finally, the 
requirement r tr <C t 9 = 1/J^/p leads to the condi- 
tion J <C J t " ~ d H 2 /y/p. In macroscopic systems the 
quantities and J" r are very large (for example at 
£ ~ lO-^m and ~ / L ~ 10a we have J£ - 10 14 
and J t ™ ~ lO 12 /^). At temperatures of intensive des- 
orption (n/a 2 Dn > 1) the quantity Jg is as large. So, 
the reduced flux density J can, formally, be taken as be- 
ing unlimitedly variable. 



III. TRANSITION IN RELAXATION DYNAMICS 

At a high enough density of the input flux, J, a 
quasiequilibrium J d — h s l s ~ J should be established 
well before the particles distribution in the bulk becomes 
uniform. In view of the fact that at p < 1 there must 
be l s < h s , this means h s /h e > 1, i.e. at sufficiently 
large J the H particles surface concentration should ini- 
tially grow to some maximum hf 1 > h e and then, as a 
result of diffusion into the bulk, relax asymtotically to 
its steady state value, h e , from above. We thus conclude 
that at p < 1 there ought to exist a critical flux density 
J c (p), above which the relaxation character is qualita- 
tively changed. In this section we give a linear analysis 
of the long-time relaxation dynamics, in terms of which 
we find the relaxation rate us(J,p) and the exact value of 

We introduce the new variables h — h e — h,l = l e — I 
and h = N e — N. Then, instead of (6) we have 

dh/dT = V 2 h, dl/dr = (l/p)W 2 !, (7a) 



Vh \ s = (l/p)V/ | s = -VJ(h. + I) + Us (7b) 

with V(/i, I) | x= o= 0- In a long-time limit h,l | T ^oo^ 
the nonlinear term may be neglected, so by satisfying 
(7) , we find the leading terms of the long-time relaxation 
in the form 

h = Ah cos(\fux)e~ UJT ', I — Al cos( v / pcUx)e~" r , (8) 

where Al = (^sin^/uj/ sin ^/puj)Ah , and the relax- 
ation rate iv(p 7 J) is defined by the least positive root of 
the equation 
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cot -Juj = — J~p~ cot y/puj 



(9) 



It follows from (8) that long-time asymptotics of sur- 
face concentrations as a function of the particles number 
h = Jq 1 hdx = Jq Idx has the form 

h s |n^o= cot y/w)n ,l s |ft^0 = {y/puj cot ^/puj)h, (10) 
whence 



\A^0- 



cot \fuj 
Jpcot^pHJ 




tan(- s /paJ). (11) 



According to (9), with growing J the relaxation rate is 
increased from to = 2\f~J in the reaction-controlled limit 
VJ< 1 

w = 2>/j(l--^-^Vj+---), \/j<l 



to a maximal 

w = w m (p)(l - b(p)/VJ ■ 



\/j> 1 



in the diffusion-controlled limit 3> 1, the maximal 
relaxation rate Lo m (p) growing herewith from 



W m (l)=7T 2 /4 



at p — > 1 to 



calculated from (9), (11), the second one demonstrates 
the time dependences of h s /h e and l s /l e at J = J c and 
J = 10 3 , numerically calculated from (6). 

Eqs. (10) and (11) carry no information on the initial 
conditions, so they, clearly, remain valid at any initial 
number of particles A//r(0) = A/z,(0) = 7V(0) and at their 
arbitrary initial distribution. In this case, according to 
(11), if the initial distribution is uniform, and n > 0, i.e. 
the initial particles number is less than the steady state 
one, then at J > J c the h s value passes a maximum and 
relaxes to the steady state level from above, whereas if 
the initial distribution is uniform, and n < 0, i.e. the 
initial particles number is greater than the steady state 
one, then at J > J c the h s value passes a mimimum and 
relaxes to the steady state level from below. 

In the present paper the main attention will be focused 
on the behavior of the surface concentrations and the 
growth dynamics of the particles number in the diffusion- 
controlled limit J ^> J c , which is of greatest interest. As 
at lot > 1 deviations from the steady state become expo- 
nentially small then, clearly, in the diffusion-controlled 
limit the main nonlinear kinetic effects evolve in times 
r <C 1 when the diffusion of H particles proceeds, ac- 
tually, into a semi-infinite medium. In the next Section 
we shall consider the dynamics of behavior of h Sl l s and 
N at r -C 1, taking into account the reflection from the 
boundary x = 0, and then, in Section V, we shall come 
back to an analysis of the long-time relaxation and obtain 
expressions for the Ai(J,p) amplitudes. 



L0 m {0) = 4.1158... 

at p — > (b(p) changes from 6(1) = 1 to 6(0) = 
1.348...). So, in accord with (9), (10), we conclude that 
at p < 1 there exists a critical relaxation rate u> c — 7r 2 /4, 
which is reached at a critical flux density 
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Jc = '—t'AT?{TT^/p/2), 



(12) 



above which the quantity h s reverses its sign (H ► — ), 

i.e. the character of the surface concentration relaxation 
of H particles qualitatively changes: at J < J c the h s 
value grows to h e monotonously, whereas at J > J c the 
value of h s reaches initially a maximum and then relaxes 
to h e from above (Fig. 2). According to (11), as J grows, 
the ratio r s = h s /l s |„_>o is changed from 1 at ^/J/J c <C 



r a = l--VJ(l-p) + .. 



\P7Tc^ i 



to -1 at ^JJ/Jc > 1 



r s = -1 + 




tan y/pu> m , y/J/Jc > 1 



IV. DIFFUSION OF H PARTICLES INTO A 
SEMI-INFINITE MEDIUM (t«1). 



By applying the Laplace transform f(s) = Cf{r) = 
Jo°° e_ST /( T )^ T to Eqs. (6), one can easily obtain a for- 
mal solution of problem (6) in the form 

w x ? cosh(xyi) - j cosh(x^/sp) 

h(x,s)=hs — l(x,s) = l s - (13) 

cosh Jap 



cosh-yi 

with the boundary conditions (6b), which can, for con- 
venience, be represented as 



N = —]= tanh ^fs = — ^= tanh Jsp 
\fs ' J^P 



and 



N = ^-s- 1 C(h s l s ). 



(14) 



(15) 



The system of equations (14) and (15) comletely de- 
fines the behavior of N(r),h s (T), and I s (t) which, in 
turn, define the evolution of spatial particles distribu- 
tion via Eqs. (13). In the limit of our interest here r«l 
the first of Eqs. (14) is reduced to 



as is illustrated for p = 0.1 in Figs. 3 and 4, the first 
of which demonstrates the dependences of r s on J/J Cl 



N 



2e 



+ •••). 



(16) 
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where in expansion of tanh y/s in power-series of e~ 2 ^ 
the leading term 0(e _1 ^ T ), which takes into account the 
contribution of reflection from the boundary x = 0, is 
retained. The second characteristic time tl = p, at which 
the diffusive length L of particles becomes equal to the 
system's size, separates two temporal regions, T<p and 
t » p, within which the growth dynamics of the particles 
number is qualitatively different. Within r <C p the L 
particles duffusion proceeds, actually, into a semi-infinite 
medium, and from (14) we have 



N= -^(1-26-^ + ...). 



(17) 



In the opposite limit, t>ji, the L particles distribu- 
tion becomes practically uniform, and from (14) it follows 



N = l s ( 1 - -sp + ■ ■ 



(18) 



where in expansion of tanh ^/sp in power-series of y/~sp 
the leading term 0(p/r) is retained, which in the case of 
the power growth l s oc r™ leads to the law 



N = l s 



(19) 



In what follows, we shall consider the dynamics of 
crossover from the adsorption-controlled regime to the 
diffusion-controlled one first in the limit of t <C p, when 
H and L particles diffuse into a semi-infinite medium, 
then in the limit p <C t <C 1, when H particles diffuse 
into a semi-infinite medium at a uniform L particles dis- 
tribution. 



A. Diffusion of H and L particles into a semi-infinite 
medium (r <C p). 

On neglecting the reflection from the boundary x = 0, 
we have from (16) and (17) 



(20) 



The system of equations (15) and (20) is reduced to 
the nonlinear integral equation that cannot be solved in 
the general form. Our aim will be, starting with (20), 
to obtain asymptotics for the solution of (15)-(17) in the 
limits of J D | s ~ J » J d (adsorption-controlled regime) 
and J d ~ J » J D \ s (diffusion-controlled regime). 



whence it follows that the adsorption-controlled 
asymptotics has the form 



h s = —=Jy/r(l + q h ), 
V 71 " 



V 7r 



N = Jr(l + g N ), 



(21) 



where ft< 1. Thus, the condition of smallness of des- 
orption flux is J d — h s l s — (4/7r)J 2 ^/pr <C J, whence 
it follows that it takes place at r <C r q = 1/J^/p. 
Substituting the expression for J d into (15), we obtain 
N = Jr(l - |(r/r g ) + •••). Substituting further this 
expression into (16) and (17), with an accuracy to the 
leading terms, we finally obtain 



Qh = ~7T ( — J + 2Vnierfc—= + 

37T \T q J y/T 



(22) 



0N 
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+ 



where the function ierfc(rf) = e v / 'y/n 
at t] 1 has the asymptotic form e~ n (1 
-)/2y/W- 



- rjerf(r)) 
3/2T7 2 + 



2. Diffusion- controlled regime (r q <JC r <§C p). 

In this limit, by neglecting the contribution of the tran- 
sient region, i.e. by assuming r q — > 0, from (20) and (15) 
we have 

h s = l s /\/p, h s l s = J, 

whence we find h s = h e p~ 1 / 4 ,l s = l e p x l 4 and, 
subsequently, according to (20), we obtain N = 
2 v /jA' : p _1/4 \A : and, hence, J D \„= J(^) 1/2 /y/n. Thus, 
the diffusion-controlled asymptotics has the form 

h s = h e p- 1 '\l + \ H ), 



1. Adsorption- controlled regime (r -Cp, r 9/ ) 



l s =l eP 1 '\l + \ L ) 



(23) 



By assuming N = Jt, i.e. by neglecting the contribu- 
tion of desorption, from (20) we obtain 

h s = l,/y/p = (2/V^)JV7, 



N = -^N e p-^ 4 V^(l + Xn), 

v 71- 
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where A; < 1. Substituting (23) into Eqs. (15), (16), 
and (17), we find that, asymptotically, Xh and Al are 
connected by relationships 



A/, - A- - : -2 ( ; i JcJZ - erfc-^= 



1 / T \ !/ 2 

\/tt V r / 



whence, with taking account of (16) and (17), we fi- 
nally obtain 

A tf = -7T^(— -erfcJ-+erfc-=, 

Zy/TT \ T ) \ T y/T 



1 /T\l/2 /77 i 

■2^(7) +-/V"- e ^ c -A' (24) 



A 



X (7) V2 " ^ ( 2er/c /f + ierfc ^) 



where comlementar error function erfc(rj) = l — erf(rj) 
at 77 <C 1 has the asymptotics e~ v (1 — l/2i] 2 + ■ ■ ^/y/nrj. 



I s = Jt(1 + A l ), 



N = Jt(1 + Ajv). 



(27) 



According to (26), the condition of smallness of the 
desorption flux takes now the form J d = h s l s ~ 
(2/y/Tr)J 2 T 3 / 2 <C J, whence it follows that it is real- 
ized at times t<t, = J~ 2 / 3 . So, under the conditions 
of uniform distribution of L particles, there appears a 
new characteristic time scale r* = J~ 2 / 3 in the vicinity 
of which the crossover from the adsorption-controlled to 
the diffusion-controlled regime occurs. Using the exact 
series expansion in particular points 



coth y/sp 



sp 



y/sp ^ sp + (irn) 2 ' 



we find from (14) 



(- 00 
l„=N{l + -ap-2J2 

{ 6 n=l 



(sp/irn) 2 



(28) 



sp + (7m) 2 j ' 
and, therefore, in the limit r/r* — > we have exactly 



=Jr 1 



(?) 



1 
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n=l 



n 2 e 



B. Diffusion of H particles into a semi-infinite 
medium at uniform distribution of L particles 

(p < T < 1). 

By neglecting the H particles reflection from the 
boundary x = and nonuniformity of the distribution 
of L particles, from (16) and (18) we have 



N = l s = h s /sfs. 



(25) 



whence it is seen that the addition to (19) becomes 
negligibly small 0(e _7r ) already at r <~ p. Substituting 
the expression for J d into (15), we obtain N = Jr(l — 
j^(-j-) 3 ^ 2 + •••)■ Substituting then this expression into 
(16) and (28), respectively, we finally find 



A 



H 



8 V^* 



3/2 



2y/nierfc 



1 



Evidently, the character of the crossover onto regime 
(25) depends on the t q /tl relation. In the region of r q <C 
tt, regime (25) is realized after a quasiequilibrium h s l s ~ 
J has been reached, whereas in the opposite limit, tl <C 
T q , the crossover to regime (25) occurs directly at the 
initial stage when the desorption can yet be neglected. 



A H©-5^(if 2+ 



4 / T 



By/lT \T* 



3/2 



+ 



(29) 



1. Adsorption- controlled regime fp « t « r* ). 

By assuming N — Jt, i.e. by neglecting the contribu- 
tion of desorption, from (25) we have 



N = l s = Jt, h s = {2/y/^)Jy/^, 



(26) 



whence it follows that the adsorption-controlled 
asymptotics has the form 



h s = ^Jy/T(l + A H ), 



2. Diffusion- controlled regime (p, t* <C t <C 1) 

By neglecting the contribution of the transient region, 
i.e. by taking p, t* — > and assuming that /i s changes by 
the power law h s <~ r™, from (25) and (15) we have 



N = l s = j(n)h s y/T, h s l s = J, 



(30) 



where j(n) = T(n + 1)/T(n + 3/2) and T(n) is the 
gamma function. From (30) it immediately follows n = 



6 



— 1/4. Thus, the exact asymptotic solution of (15), (25) 
in the limit p — > 0, J — > oo is 



N = l s = Z^r 1 / 4 , h s = hrf-W 



-1/4 



(31) 



where (3 = [r(3/4)/r(5/4)] 1 / 2 = 1.1627366.... With 
account taken of the contribution of the transient region 
and reflection of H particles from the boundary x = 0, 
we, finally, obtain the diffusion-controlled asymptotics in 
the form 



l s =l e /3T 1 / 4 (l + ( j )L ), 
N = N e f3T 1/4 {l + (f> N ). 



(32) 



Substituting (32) into (15) and (18) and taking <fo <C 1 
we find that, asymptotically, 4>h,4>l, and <f>N are con- 
nected by the relationships 



bH + 4>L = - 1 {-) + 



and 



:J N 



(33) 



(34) 



Substituting then (32) into Eq.(16), we come to the 
equation 



e -v", (35) 



which, on neglecting the contribution of H particles 
reflection from the boundary x = 0, takes the form 

pCtfNT 1 ' 4 ) = s-Wcfa/T 1 ' 4 ). (36) 

Assuming that <pi decay asymptotically in time by the 
power-law 



= Ci/T Ui 



(37) 



where C»(p, r*) — > as p, r» — > 0, from (33), (34) and 
(36) it can easily be shown that the exponents V{ > 3/4. 
Indeed, let uh < 3/4. Then, from (36) it follows that 
ul = vh, and the coefficients Cl and Ch have like signs, 
which contradicts the condition (33). Thus, i>$ > 3/4 
which suggests that in contrast to the diffusion-controlled 
asymptotics (23), the sign and value of the Cl and Ch 
are governed by all the prehistory of transition to the 
asymptotics (31) and can be calculated only numerically. 
In the next sections we shall discuss the behavior of <j)L 
and <Ph in more detail and demonstrate that it is quite 
nontrivial. We now consider the contribution into <pi due 
to the reflection of H particles from the boundary x = 0. 
According to (33) and (34) , in the limit p, t* — > we have 



<j>H = -<f>L = ~<f>N = ^ (0) - 
Substituting <fi^ into Eq.(35), we obtain 

pt&W) = ^f^e- 2 ^ - S - 1 /2£(0(o) /r V4 ) . (38) 

It is easy to check that the inverse Laplace transform 
of the solution of Eq.(38) ought to have the form 



^ = « T 3/4 e -l/r (1+s)i 



(39) 



where 



9 = c 1/2 t 1/2 + cir + c 3/2 t 3/2 + 



Substituting Eq.(39) into (38) and equlizing the coef- 
ficients at the same s powers, we finally find 



1 



15 1 



1-ci 



C 1/ 2--^,C 1 --- + —,C 3/2 ^ 



C. Two characteristic paths of the crossover to the 
diffusion-controlled regime with uniform L particle 
distribution. 

According to subsections A) and B), depending on the 
relation of three characteristic times tl = P,T q , and r«, 
the crossover to the diffusion-controlled regime with uni- 
form L particles distribution proceeds by one of the two 
qualitatively different scenarios. One of these is realized 
in the limit T q «p« 1, when at first, a quasiequilib- 
rium between the input and dcsorption fluxes is estab- 
lished in the system (crossover Al — > A2 at r ~ r 9 ), 
following which the L particles distribution becoms uni- 
form (crossover A2 — > B2 at r ~ p). In this case, the 
surface concentrations h Sl l s and particle number N go, 
respectively, through the following chains of power-law 
asymptotics: 



h a 



,1/2 



-1/4 



-1/2 



-1/4 



(40) 



N : 



rV2 



-1/4 



The second scenario is realized in the limit p « t* « 1, 
when at first the L particles distribution becomes uni- 
form (crossover Al — > Bl at r ~ p), whereupon a 
quasiequilibrium is established between the input and 
desorption fluxes (crossover Bl — > B2 at r ~ r*). In 
this case, h s ,l s , and N go, respectively, through the fol- 
lowing chains of power-law asymptotics: 



-1/2 



-1/4 



7 



.1/2 



N 



-1/4 



,1/4 



(41) 



It can easily be shown from Eqs.(13) that at each of 
the given power-law asymptotics, h s oc r",r <C 1, the 
spatial distribution of iJ-particles changes, on neglecting 
the contribution of the transient regions, by the law 



where 



h(x,r) = h s T n (tt^t) 



( ir^ierfciO, 



(42) 



•MO = { 



n = 1/2 



erfc(£), 

^zi) e -« 2 / 2 p_ 1/2( v2e), 







-1/4. 



Here T> v (z) is the function of parabolic cylinder which 
at z » 1 has the asymptotics T> v (z) = e~ z / 4 z v [l — v(y — 
l)/2z 2 + ...]. The spatial distribution of L particles at 
asymptotics l s cx t™,t -C p with n = 1/2,0 has exactly 
the same form, except for the substitution t — > r/p. 

By comparing the paired ratios of the characteristic 
times TL,r q , and r», wc find 



tl 



3/2 



K 



3/2 



J, 



(43) 



whence it is seen that the character of the system's 
evolution is governed by the value of the parameter K = 
p 3 / 2 J: at any p and J relation in the limit /C > 1 the 
system evolves in accord with the chains of asymptotics 
(40), whereas in the opposite limit, K, <C 1, the system 
evolves in accord with the chains of asymptotics (41). On 
neglecting the reflection of H particles from the boundary 
x = 0, one can easily see that in the limit JC — > oo(r q /p — > 
0) the only characteristic time scale is defined by the 
quantity p, therefore, as suggested by (23) and (32), the 
behavior of h s and l s should have the scaling form 



h s /h e = p- L/i w H (T/p), l s /l e = P 1 '*w l (t/ P ), (44) 

where the scaling functions Wi(() have asymptotics 
Wh,l{0 = 1 T er/c(l/VC) at C « 1 and w h ,l(Q - 
(/?C 1/4 ) Tl at C > 1. Indeed, substituting (44) into (14) 
and (15) and going to the Lasplace transform with re- 
spect to the reduced time Q = r/p —f s, in the limit 
JC — > oo wc find 

whWl = 1) w H = w L tanh-^/s, 

whence it immediately follows Wi = Wi(r/p) and, 
hence, tfii = 4>i{ T /p)- I n the opposite limit, /C — > 
0(p/r» — > 0), the only characteristic time scale is defind 
by the quantity r*, therefore, in accord with (27) and 



(32), the behavior of /i s and / s should have the scaling 
form 

h s /h e - J 1/6 v h (t/t,), l s /l e = J- 1/6 v l (t/t,), (45) 

where the scaling functions i>i(£) have asymptotics 
Wff(C) ^ 2^C/^,v L (C) ~ C at C < 1 and VffiL ~ 
(/3C 1/4 ) Tl at C > 1. Indeed, substituting (45) into (14) 
and (15) and going to the Laplace transform with respect 
to the reduced time £ = r/r* — > s, in the limit /C — > wc 
find 

£{v H v L ) = s _1 (l - s 2 «l), v h = v l ^, 

whence it immediately follows = Uj(r/r*) and, 
hence, ^ = ^(r/r*). 

From (45) it follows that in the limit /C — > the 
crossover to the diffusion-controlled regime is character- 
ized by a comparatively sharp maximum of h s /h e , the 
height of which depends solely on J and changes with 
growing J by the law 



hf/h e = bnJ 1 '*. 



(46) 



In the opposite limit tC — > oo, according to (23), (44), 
a maximum of h 8 /h e is degenerated to an extended 
plateau, of which the height depends on p alone 



K/h e = P 



- „-V4 



(47) 



and the extension, defined as the ratio of the times, 
which bound the region A# < e <C 1, grows with J as 
oc/Ce 2 /ln(l/e). 

In the next subsections we give the results of a detailed 
numerical study of the dynamics of the system in ques- 
tion, enabling one to gain a complete picture of evolu- 
tion of h s ,l s , and TV "trajectories" between two limiting 
regimes (40) and (41). 



D. Numerical calculations 

The numerical integration of equations (6) was per- 
formed by means of the implicit discretization scheme of 
increased accuracy with an additional "fictitious" node 
at the surface. The scheme allowed performing the cal- 
culations in the system with strong difference in species 
diffusivities with an accuracy down to 10^ 3 % (see below). 
The space and time steps were changed within the ranges 
5x ~ 3x 10~ 4 ^3x 10~ 6 and St ~ 10" 4 10" 11 , respec- 
tively, with the number of time steps being 10 5 -^10 6 . The 
evolution dynamics of /i s (r), I s (t), and N(t) was studied 
in the ranges of p = 10~ 10 1 and J = 10" 1 + 10 12 . 



1. Kinetic diagrams p — J — r 

Fig. 5(a) illustrates 
the time dependences h s (T)/h e ,l s (T)/l e , and N(r)/N e , 



8 



calculated numerically at a fixed flux density J = 10 8 
for p values ranging from f to 10 -9 (the arrows show the 
directions of shift of the corresponding trajectories at a 
change of p from p = 1 (bold curves) to p = 10~ 9 (bold 
curves)). One can clearly see the whole set of power-law 
asymptotic regions (40), (41) of h s ,l s , and N trajectories 
and the character of their evolution as p (and, hence, JC) 
are decreased from the values corresponding to the lim- 
iting regime (40) (p > 5 x 10~ 6 ,/C > 1) to those of (41) 
(p <C 5 x 10~ 6 ,/C <C 1) where the dependences of the 
trajectories on p dissapears and they asymptotically ap- 
proach those shown in bold lines. In accord with (46) and 
(47), as p decreases, the plateau on the curves h s {r)/h e 
is transformed to a comparatively sharp maximum, the 
height of which depends only on J. At all p <C 1 the 
h s /h e ,l s /l e , and N/N e trajectories come to the univer- 
sal asymptotics (/^t 1 / 4 )^ 1 (31). Herewith, as seen from 
Fig. 5a and will be demonstrated in detail in what follows, 
the h s (r)/h e trajectories always come to the asymptotics 
(/St 1 / 4 ) -1 only from above, whereas the / S (r)/Z e trajec- 
tories always come to the asymptotics /3-r 1 / 4 only from 
below. 

By defining the boundaries of the asymptotic regions 
so that within their confines the condition 

max{|(5 |, \d5 a /d\nT\} < 0.01 (48) 

shall hold (here S a = {ft, \i, Aj, fa, crj}, Ui describes 
the exponential relaxation to the steady state according 
to Eqs.(69)), we have carried out an extensive numerical 
study of the positions of the corresponding boundaries 
for J = 10 8 and p values, ranging from 1 to 10~ 10 (we 
have calculated and analyzed trajectories for 100 values 
of p, 10 for each order). Figs. 5b and c show the kinetic 
p — t diagrams of the regions of the power-law asymp- 
totics and steady state (shaded by light gray) for the 
h s {T)/h e and l s (r)/l e trajectories, respectively. In Fig. 5c 
by dark gray is distinguished the /Jt 1 / 4 asymptotic region 
for the N(r)/N e trajectories of the growth of the particles 
number. From Figs. 5b, c it is seen that the dashed lines of 
the characteristic times Tq, , tl — p and th — 1 divide 
the p — t plane into several segments, confining the re- 
gions of the power-law asymptotics r™ with n — 1/2, 0, 1 
and Tl/4- The boundaries of these regions go in par- 
allel with the lines of the corresponding characteristic 
times, in accord with Eqs. (22), (24), (29), (39), (44) 
and (45). The t q , tl , and r* lines intersect in the point 
JC = 1 (r g = t l = t* = 4.64158 x 10~ 6 ), shown by filled 
circle, in accord with Eq.(43). As the n = 1/2,0,1 re- 
gion boundary positions are described by Eqs. (22), (24) 
and (26), we shall focus mainly on the n = =Fl/4 region 
boundary positions. 

From Figs. 5b, c it is seen that, in accord with (44), 
at p » 5 x 10~ 6 (/C S> 1) the left-hand boundaries of 
the n = =Fl/4 regions go in parallel with the tl line 
down to the point of intersection of left- and right-hand 
boundaries of n = regions (p«2x 10~ 3 , JC w 5 x 10 3 ) 
where the plateau on the h s and l s trajectories disap- 
pears. Above this point the h s /h e and l s /l e approarch 



the (/St 1 ' 4 )^ 1 asymptotics simultaneously and, which is 
important, this occurs at r w p, i.e. long before the L 
particles distribution in the bulk becomes uniform (left 
boundary of the n = 1/4 region shaded by dark gray 
in Fig. 5c). Such synhchronization of the trajectories di- 
rectly follows from Eqs. (33) and (44), according to which 
in the scaling limit JC — > oo (r*/f> — > 0) whwl = 1 and 
4>h = —<Pl- In the opposite limit p«5x 10~ 6 (/C <C 1) 
the left boundaries of n = Tl/4 regions become par- 
allel to the t* line in the vicinity of the intersection 
point of n = 1/2 boundary (Fig. 5b) with the tl line 
(p w 3 x 10~ 7 ,/C w 2 x 10~ 2 ), in accord with Eq.(45). 
Below this point (i) the behavior of h s (r) obeys the scal- 
ing law (45), i.e. it ceases to be dependent on p and 
(ii) the boundaries of n = 1/4 regions of growth of 
I s (t) and N(t) merge, i.e. by the instant of the tran- 
sition to asymptotics /3-r 1 / 4 the distribution of L parti- 
cles becomes uniform. Essentially, that in this region 
of p and JC the behavior of h s and l s becomes strongly 
asymmetric: the transition of the h s /h e trajectories to 
(/3-r 1 / 4 )" 1 asymptotics takes place well before the l s /l e 
trajectories reach the fir 1 ' 4 asymptotics. Below the point 
p w 2 x 10~ 9 ,/C « 10 -5 , where on the kinetic diagram 
of Fig. 5c there appears the n = 1 region, the merging of 
the I s (t) and N(t) trajectories takes place already on the 
stage of the adsorption-controlled regime, beyond which 
the l s (r) trajectory ceases to be p- independent. 

With account taken of the fact that the boundaries of 
the asymptotic regions are strongly related to the char- 
acteristic times TL,T q , and r*, one can easily get from 
Fig. 5b, c the idea of how the kinetic diagrams evolve with 
changing density of the input flux J. Indeed, as J grows, 
the point JC = 1 should shift along the tl line, as shown 
by arrows in Fig. 5b, c, entailing the T q and t* lines. The 
T q and r* lines, shifting in parallel to themselves, should 
in turn entail the lines of the corresponding boundaries 
of the asymptotic regions, so that with growing J the 
n = 0,=Fl/4 regions should expand in a self-similar man- 
ner, and the n = 1/2, 1 regions contract in a self-similar 
manner. From the conditions r_i/4 ~ r r and n/4 ~ r r , 
where t_!/ 4 = 5.855-r* and T\u — 93.88-r* are the po- 
sitions of the left boundaries of regions n = —1/4 and 
n=l/4at/C^0 and r r w 0.2 is the position of their 
right boundary associated with the influence of the H 
particles reflection from the boundary x — 0, one can 
conclude that the starting flux densities for the appear- 
ance of t -1 / 4 and r 1 / 4 asymptotics are J-1/4, w 10 2 and 
J1/4 s=s 10 4 , respectively. From Fig. 5b, c it follows that 
the vertices of n — Tl/4 regions for the surface concen- 
trations and the number of particles are in the points 
p s rts 0.1 and pfq rts 2 x 10~ 2 , respectively. Thus, the 
J — p limits of (/3r 1 / 4 ) Tl asymptotics are defined by the 
following conditions: 

h s :p<p s K 0.1, J > J_ 1/4 w 10 2 , 
l s :p<p s K 0.1, J > J 1/4 w 10 4 , 
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N : p < p N w 2 x 10~ 2 , J > Ji/4 « 10 4 . 



2. Maximum of h s 

Figs. 6a and b illustrates the behavior of the maximum 
of the surface H particles concentration, h^f /h e ,(&) and 
the time for which this maximum is reached, t M>(b) as 
a function of growing J. The curves are calculated nu- 
merically at fixed p = 10 _1 ,10- 2 ,10- 3 , 10" 4 , 10" 5 and 
10~ 6 . It is seen that, in accord with (45), (46), at small 
JC <C 1 the calculated curves approach the p-indepcndcnt 
power-law asymptotics (shown in the dashed lines) 

hf /h e = 0.7221J 1 / 6 , t m = 1.125J" 2 / 3 . (49) 

As J grows, the dependences of /h e and tm on J 
deviate from the p-independent asymptotics (49) the ear- 
lier the large is p, reaching at JC 3> 1 the J-indepcndcnt 
plateau (47) (as seen for the curves for p = 10 _1 and 
p = 10~ 2 ). Interestingly that at p <C 1 the curves ap- 
proach the power-law asymptotics (49) already at com- 
paratively small flux densities J w 50, exceeding J c by 
no more than an order of magnitude. 

3. Behavior of 4>h and 4>l 

We shall now consider the regularities of the behavior 
of transient terms 4>h{t) and 4>l(t), which characrerize 
the kinetics of the transition of surface concentrations 
to the universal asymptotics (/Jt 1 / 4 )^. Fig. 7 shows the 
dependences 4>h(t){&) and 0L(r)(b), calculated numeri- 
cally at p = 10~ 4 for J values, ranging from J = 10 4 to 
J = 10 12 (from JC = 10~ 2 to JC = 10 6 , respectively). It 
is seen that (i) the <pjj value at any JC first crosses zero, 
changing the sign from - to + , then reaches a maximum 
and asymptotically approaches zero only from above; (ii) 
the <pL value at any JC asymptotically approaches zero 
only from below, in this case, with JC < 1 it occurs 
monotonously whereas with C > 1 the <f>L value first 
crosses zero, changing the sign from + to - , then reaches 
in modulus a maximum, and only after this begins to 
approach zero; (iii) with a growth in J the 4>h maxi- 
mum is shifted left and its amplitude grows, whereas the 
\(j>L\ maximum is shifted right and its amplitude drops, 
so that at JC > 10 6 the behavior of 4>h(t) and 4>l(t) be- 
comes comletcly "symmetric", in agreeement with (33), 
(44) 

Mr) = \Mr)\- 

From Fig. 8 where are presented the sections of the de- 
pendences 4>h(t) > 0(a) and </>l(t) < 0(b) replotted in 
double logarithmic coordinates, we find that at J » Ji/ 4 
the <j>H and 4>l values decrease at a sufficient distance 



from the <pjj maximum by the power law (37) with the 
exponents vjj = Vr, = 3/4 

4>i = C,/t 3 / 4 + ■■■ (50) 

up to r rs 0.1 -j- 0.2 where a rapid growth of \<j>i \ begins 
due to the reflection of H particles from the boundary 
x = 0. At the growth stage the \cj>i\ dependences are 
seen to go exactly onto the bold (j>^ curve, calculated 
from Eq.(39). It is important to note that from Fig. 8 it 
directly follows that the error of numerical calculations 
does not exceed ss 10~ 3 %. According to Eq.(44) in the 
scaling limit JC — > oo from (50) and (33) it follows 

^mr(p/r) 3/4 + '", 

where the coefficients — —m^ are independent of 
J and p. From the data of Fig. 8 we find at JC > 10 6 

= -m^ = 0.023 (51) 

and obtain for the maximum \(f>i\ m and the time of its 
attaining r™ 

0£ = \<j, L \ m = 0.00798, = 2.51p. (52) 

In Fig. 9 are plotted in double logarithmic coordi- 
nates the sections of the dependences 4>h(t) > 0(a) and 
4>l(t) < 0(b), calculated numerically at a fixed flux den- 
sity J = 10 6 for the p values, ranging from p = 10~ 2 to 
p = 10~ 6 (from JC = 10 3 to JC = 10~ 3 , respectively). The 
both maxima, <j)jj and |</>i,|' M \ are seen to shift towards 
smaller r with decreasing p, herewith the 4>h maximum 
amplitude decreases, whereas the \4>l\ maximum ampli- 
tude increases, and at JC < 1 it disappears. At p <C p s 
the 4>h{t) and 4>l(t) dependences decay at a distance 
from (jff^ by the power law (50), converging to the p inde- 
pendent trajectories at p w 10~ 6 (/C ss 10~ 3 ). According 
to Eq.(45) in the scaling limit JC — > from (50) and (33) 
it follows 

^m°(r,/r) 3 / 4 + ..., 

where the p- and J-independent coefficients m° should 
satisfy the condition 

m° H + m° L = -13/ A = -0.29068.... (53) 

From the data of Fig. 9, in accord with (53), we find at 

JC < 10- 3 

m° H = 0.014, m° L = -0.305 (54) 

and obtain for the <fff!j maximum and the time of its 
attaining 

<j>% = 0.00108, t% = 17.3r,. (55) 
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E. Scaling and universality 

So far the main attention has been focused on the lim- 
iting regimes K < 1 and > 1 and on specific features 
of the crossover between these regimes at a change in J 
for p =const. or at a change in p for J =const. We shall 
now show that at r <C 1, when the H particles reflection 
from the boundary x = can be neglected, the behavior 
of h 8 ,l s , and N is described by the scaling laws 

h s =h e p- 1 / i W H (K,T), 
I s =l e p 1 / 4 W L (1C,T), (56) 

N = Nep^WwilCT), 

where T = r/p and the scaling functions Wi(JC,T) 
depend on the only parameter tC — p 3 / 2 J. Indeed, going 
to the Laplace transform with respect to the reduced time 
T — > s, on substituting (56) into Eqs.(14) and (15) and 
neglecting the terms, related to the H particles reflection 
from the boundary x = 0, we easily find 

Wn = -^ = ^ tanh v/i, (57) 
V s V s 



B) Diffusion of H particles into a semi-infinite medium 
at uniform distribution of L particles (1 <C T <C 1/p)- 

1. Adsorption-controlled regime (1 <C T <C 1//C 2//3 ). 

W H = ^=VKT (l h^^ T V2 + . . )j 
W l =tVk(i + ^- ^KT 3 ' 2 + •••), (61) 

w N = tVjc (\ - ^=/cr 3 / 2 + ■ ■ ^ . 

2. Diffusion- controlled regime (T > 1, 1//C 2 / 3 ). 
W h ,l = (/JT 1 / 4 )^! + m H ^ L (IC)T- 3 / 4 + •••), (62) 

W N = (3T^ (l - ^ + m L (/C)T- 3 / 4 + -..). 

According to Eqs.(56), at arbitrary changes in J and 
p, with /C =const. retained, the scaling functions are the 
universal functions of the reduced time, Wi(T), and the 
h s /h e , l s /l e , and N/N e trajectories are shifted by a factor 



fiC 

W N = *—(1-sC(WhW l )), 



(58) 



whence it follows W, = Wi(fC,T). In agreement with 
the results of Sections IVA and B, the asymptotic behav- 
ior of the scaling functions Wi(K,T) on the power-law 
portions of the h s ,l s , and N trajectories have the follow- 
ing form: 

A) Diffusion of H and L particles into a semi-infinite 
medium (T < 1). 

1. Adsorption- controlled regime (T <C 1, 1/JC). 



W H = ^VtCT [ 1 



_8 

3tt' 



ACT- 



VK L = -^VO ; ^l-^/Cr + Te- 1 / T + --^, (59) 



2. Diffusion- controlled regime (1/JC <C T <C 1). 
Wh.l = 1 



— / =F er/c— j= + 
2\A^T VT 



W N = 2 VT> (l - \ ^ V^ierfc-L + 



(60) 



= (J/Kf 1 ^. 

So, the boundaries of the power-law asymptotic re- 
gions, the Wh and \<j>i\ maxima, the times at which they 
are attained, Tm and X™, and the rrii coefficients are un- 
ambiguous functions of the parameter tC. Their behavior 
will be considered in what follows. 



1. Behavior ofWi(K.,T) and kinetic diagrams K, — T 

Figs. 10, 11, and 
12 show the plots W H = p 1 ' i h s /h e ,W L = p-^l./U, 
and W N = p-V^JV/JVe vs T = r/p in the K range from 
K = 10~ 11//2 to K, = 10 13 / 2 (upper panels), and the cor- 
responding K — T diagrams of the power-law asymptotic 
regions (lower panels), obtained by replotting the data 
of Fig. 5 (J = 10 8 ) for r < r r 0.2 in the corresponding 
scaling coordinates. It is seen that with the growing JC 
from JC <C 1, the Wi(tC,T) trajectories initially shift in a 
self-similar manner in accord with the scaling (45) 



W l (JC,T) ^ JC ±1/e v l (JC 2/3 T), AC-»0, 



(63) 



(note that for Wj,(/C, T) this scaling takes place only 
at T > 1 whereas at T < /C~ 2/3 , according to (28), the 
scaling Wl(JC,T) = ^/K,u(T) takes place) and then at 
/C > 1 they shift in a self-similar manner at "tails" of 
T <C 1, in accord with the scaling 
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W H , L {JC,T) = k s {JCT), 



W N (JC,T)=JC- 1 / 2 A N (1CT), 



(64) 



T M (/C) 



In 16/C ' 



converging at /CT — > oo to the universal (/C indepen- 
dent) trajectories (44) 



Wi(!C, T) = Wi(T), JCT^oo. 



(65) 



Fig. 11 clearly demonstrates the evolution of the course 
of the Wl(JC, T) trajectories as a result of competition of 
two opposite tendencies: (i) growth acceleration of l s , as- 
sociated with the reflection of L particles from the bound- 
ary x = and (ii) growth deceleration of l s , associated 
with the establishment of quasiequilibrium J d ~ J, the 
former being dominant at any JC in the vicinity of T w 1, 
the latter in the vicinity of T w T* = 1//C 2 / 3 at JC < 1 
and in the vicinity of T w = 1//C at /C > 1, shifting 
with the growing JC towards T — > 0. 

The general evolution picture of the power-law asymp- 
totic regions with the growing JC between the limiting 
regimes (40) and (41) as well as of the JC and T regions, 
in the confines of which the crossovers between power-law 
asymptotics fit the scaling laws (63), (64), and (65), are 
demonstrated by the kinetic JC — T diagrams on the lower 
panels of Figs. 10,11, and 12. According to Eqs.(56), the 
given plots yield an exhaustive description of evolution of 
the h s ,l s , and N trajectories at the stage of the H par- 
ticles diffusion into a semi-infinite medium at arbitrary 
J and p changes. In particular, at p = const, the given 
plots factually describe the evolution of the h s ,l s , and 
N trajectories with the growing density of the flux J up 
to T = T r , where T r w 0.2/p is the time, begining with 
which the reflection of H particles from the boundary 
x = becomes essential. 



2. Scaling of hf 

From Eqs.(56) it follows that at p < 1, J » J c , and 
t < r r the p and J dependences of the amplitude and 
the time tm, for which h s reaches a maximum, should 
have the scaling form 



hf = Kp-'^MilC), T M = pT M (JC). 



(66) 



In the limit of small K « 1, according to (63), (46), 
and (49), the scaling M(JC) and Tm{JC) functions have 
the asymptotics 

M(JC) = 0.7221/C 1/6 , T M {K) = 1.125/CT 2/3 . (67) 

In the opposite limit of large JC ^> 1 the quantity M(JC) 
reaches its limiting value M(JC — > oo) = 1. By differenti- 
ating (60), in this limit we find 



M{fC) = 1 - 



In 16/C 
8tt/C 



1/2 



1 + 



1 



In 16/C 



(68) 



whence it follows that M(fC) approaches 1 as 1 — M oc 
^/ln/C//C, and a rapid drop of Tm{1C) tranforms at JC 3> 1 
to a slow logarithmic one. In Fig. 12 are presented the 
numerical data of Fig. 6 replotted in the scaling coordi- 
nates p^l^h 1 ^ jh e - JC and tm/p-JC, respectively. It is seen, 
that in complete agreement with (62), the points, calcu- 
lated for different p and J values, fit the scaling functions 
M(JC) and T M {JC) which at small JC < 1 and large JC » 1 
approach, respectively, asymptotics (67) and (68), shown 
in dashes. 



3. Scaling of §™ and wii 

In Figs. 13a and b are drawn, respectively, the \4>i\ m 
and T/ 71 vs JC dependences derived from the data of Figs. 8 
and 9, and the analogous data for some other p < p s and 
J > J1/4 values. The points, calculated in a broad range 
of p and J, are seen to fit, in accord with (56), the scaling 
functions (j) m (JC) and T™(/C), which in the limit of small 
JC <C 1 approach the asymptotics (55) (dashed lines) 



= 0.00108, Th = 17.3/C~ 2/3 , 

and in the limit of large K > 1 merge and approach 
the asymptotics (52) (dashed lines) 



\<pi\ m = 0.00798, 7™ = 2.51. 

According to Eqs.(33),(50), and (56), at p <C p s and 
J > J1/4 in the range of T% < T < T r the 4>i{JC,T) 
functions ought to asymptotically decay by the law 



4>i(JC,T) = m l (JC)/T 3 / i + --- 
with the coefficients m^/C), satisfying the condition 

m H (JC) + m L (JC) = 

In Fig. 13c are drawn the JC dependences of ran and 
|tol|, derived from the data of Figs. 8 and 9, and the 
analogous data for some other p«p s and J 3> J\/a- 
The points, calculated in a broad range of p and J, are 
seen to fit the scailing functions |mj(/C)|, which in the 
limit of small JC -C 1 approach the asymptotics (54), 
(dashed lines) 

\rm(JQ\ = mf V VJC, 

and in the limit of large K > 1 approach the asymp- 
totics (51) (dashed lines) 
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rri 



H 



-TO, 



0.023. 



Figs. (7), (8), (9), and (14) give a complete picture of 
the crossovers to the universal {(3t) t1 / 4 asymptotics in a 
broad range of < JC < oo, the main features of which 
can be summarized as follows: 

1) With growing JC the cff^ amplitude grows from 
the /C-independent asymptotics (55) at JC w 10 -3 to 
the /C-independent asymptotics (52) at JC > 10 5 . The 
time of its attaining, Tjp, drops up to JC w 1 accord- 
ing to the p- independent asymptotics (55), reaching the 
/-independent asymptotics (52) in the same range of 
JC > 10 5 . 

2) With growing JC at the critical point JC m = 3.0 the 
\4>l\ maximum appears, and its amplitude then drops to 
the /C-independent asymptotics (52) at JC > 10 5 . The 
time of its attaining, T™, grows and approaches the I- 
independent asymptotics (52) in the same range of JC > 
10 5 . In the second critical point JC C = 7.3 the JC > JC C 
region appears, within which the 4>l curves cross zero, 
changing the sign herewith. 

3) Surprisingly, in the narrow vicinity of the critical 
point JC m both the time Tjj and the |tol| coefficient ap- 
proach the p- independent asymptotics (55) and (54), re- 
spectively, whereas the ran coefficient approaches the I- 
independent asymptotics (51): the formers change within 
the JC < JC m range as /C~ 2 / 3 and /C -1 / 2 , respectively, the 
latter remains constant within the JC > JC m range. 



V. LONG-TIME RELAXATION DYNAMICS 

(WT > 1). 

In the light of the results of the previous Section IV, 
we return, in conclusion, to the analysis of the long-time 
relaxation dynamics of h s ,l s , and N, which, according to 
Section III can be given as 



h a = h e (l - cr H ), 



N = N e {l-o N ), 



(69) 



where <jh = h s /h e ,aL = l s /l e , and on — n/N e at 
sufficiently large lot > 1 decay by the law 

a i = A i e- UT {l + 0{e- UT )), 

and the coefficients Ai are related as 

A H = r s A L = A N \/ui cot \/lu. (70) 

By preserving in (8) the next-to-leading terms 

h s = Ah cos(^/ujx)e~ LJT + Bh cos(V2ujx)e~ 2LJT + 



l s = Al cos(y/pujx)e UT +Bl cos(y/2pu>xe 2ulT + 
from Eq.(7) we easily find 

tn = A ie -" T (I + B i( j N + ...), (71) 
where the coefficients Bi are defind from the relations 
Bh = y/pto cot y / pwcot V2lu/Q, 

Bl = y/pto cot \fuj cot \j2pui/Q, 

Bjy = y/puj COt y/pCU COt 



Q = cot v 2uj + sjp cot 



xJ2uj/J. 



By substituting now (69) into (14) and (15) and intro- 
ducing the designation V — <t h <tl, after having excluded 
op (s) and ctjv(s), we find 



cothy^— (1 — st> (s)) y/p coth y/ps + \J 
s(cothy / s + yjpcothy/ps + \fsJJ) 



(72) 



In accord with the results of Section III, the main pole 
of Eq.(72) (s = — lo) defines the relaxation rate w as 
the least positive root of (9) and yields the relaxation 
"amplitude" Al in the form 



where 
A L o 



A l = A L0 (1 + §wX>(-w)) 

4y/p/bJ cot yjpui 
smT 2 y/u + p sin -2 y/puj + 1/ \fj 



(73) 



(74) 



and t>(— u) 



J ghGl£ u1 dr. In the diffusion- 

o 

controlled limit of our concern here, J — > oo, the coeffi- 
cients Bi change from Bn = Bl = Bh — 1 = 0.1079... at 
p — > to Bi — > at p — > 1, therefore given that the next 
after the main root of (9) is an order of magnitude above, 
we conclude that down to r ~ lj^ 1 the asymptotic be- 
havior of <jl has the form ol = A^e - " 7 (1+Bl&n + ••■) — 
Alz~^ t ■ By matching this asymptotics at r = lj^ 1 
with that of (32), which for p <C 1 assumes the form 
a L = (1 - fir 1 ' 4 ^ + 4> L )) ~ 1 - Pt 1 / 4 , in the limit p -> 
we find 



Al — e(l — (3/oj 1 / 4 ) ~ 0.5 



(75) 



From (73) and (75) it follows that at p — ► the quan- 
tity V(—lj) ~ —0.1 and the difference of Al from Alo 
does not exceed 20%. Obviously, with the growing p 
this difference can only go down, therefore by assuming 
in the first approximation Al — Alo an d extrapolat- 
ing the V ~ — A 2 L0 e~ 2l * JT asymptotics down to r = 
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(D(— ui) ~ —A 2 L0 /lj), we finally come to simple expres- 
sions 



FIG. 2. Time dependences h B (T)/h e , calculated numcri- 



A L = A H /r s ~ A L0 {1 - \A 2 L0 ), 



A N = 



4(1 



I/4 2 i 



w(sin 



psin 



pui) 



(76) 



(77) 



which, according to the numerical data, differ from 
the exact values by no more than 1 — 2% for any p. 
Fig. 15 gives, as an illustration, the time dependences 
|<7ir(r)| and ctl(t), calculated numerically at J — 10 6 for 
p = 10~ 4 , 0.5 and at J = 10 8 for p = 0.9 and compared 
with the results of the calculation from the expressions 
g i = A J e~"' T , (9), (76) and (74). The excellent agreement 
between the numerical and analytical data is evident. 
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FIG. 1. Schematic illustration of the processes of dissocia- 
tion, desorption, surface migration and diffusion into the bulk 
in the system A + B — > with input of A and B particles. 



ID At and l B = \/Dst are the diffusion lengths of A's 



cally at flux densities J = 1, J c 
(bottom to top) for p — 0.1. 



7.250031, 10^,10^, and 10 4 



FIG. 3. Dependence of r s on J/ J c , calculated from 
Eq.(5) for p — 0.1. Left and right insets demonstrate the spe- 
cific features of relaxation of surface concentrations at small 
J< J c and large J J c , respectively. 



FIG. 4. Time dependences h s (r)/h e and ls(r)/l e , calcu- 



lated numerically at J 
circles) for p = 0.1. 



J c (open circles) and J = 10 (filled 



FIG. 5. (a) Time depen- 

dences h s (T)/h e ,l s (T)/l e , and N(r)/N e , calculated numeri- 
cally at p = 1, 10 _1 , 10~ 2 , 10~ 3 , 10~ 4 , 10~ 5 , 10~ 6 , 10~ 7 , 10~ 8 , 
and 10 -9 for J = 10 8 (the arrows point to the displacement di- 
rections of the corresponding trajectories at a change in p from 
p = 1 to p = 10~ 9 ). The trajectories for p = 1 and p = 10 -9 
are given in bold lines; (b), (c) Kinetic p — r diagrams of the 
regions of power-law asymptotics and of steady state (shaded 
by light gray) for /i s (r)//i e (b) and l s (T)/l e (c) trajectories cal- 
culated numerically at J = 10 8 according to condition (48). 
The characteristic times th = 1,tl = P,T q = 1/J^/p, and 
t„ = 1/J 2//3 are shown in dashed lines. The intersection 
point of TL,r q , and r, (r q = tx = r* = 4.64158 x 10~ 6 ) is 
marked off by filled circle. The arrow shows the direction of 
displacement of this point along the tl line with the growing 
flux density J. In Fig.5c the /3r 1,/4 asymptotic region for the 
N(r)/N e trajectories is distinquished by dark gray. 



FIG. 6. Dependences of hf /h e (a) and tm (b) on J, cal- 



culated numerically at p = 10 , 10 , lO" 3 , 10 , 10" 



and 



and B's particles, respectively. 



10 (bottom to top and top to bottom in (a) and (b), re- 
spectively). Dashed lines show the asymptotics (49). 



FIG. 7. Time de- 

pendences 4>h{t) (a) and (/>l(t) (b), calculated numerically 
at p = 10~ 4 for J = 10 4 , 10 5 , 10 6 , 10 7 , 10 8 , 10 9 , 10 10 , 10 11 , and 
10 12 (right to left). 



FIG. 8. Data of Fig.7 for sections (j> H > (a) and <f) L < 
(b) replotting in double logarithmic scale. Dashed lines show 
the slope of —3/4. Bold lines show the cj>° curve, calculated 
from Eq. (39). 



FIG. 9. Time dependences 4>h(t) (a) and \4>l(t)\ (b) for 
sections 4>h > and 4>l < 0, respectively, calculated numer- 
ically at J = 10 6 for p = 10~ 2n/3 with n = 3, 4, 5, 6, 7, 8, and 
9(right to left). Dashed lines show the slope of —3/4. Bold 
lines show the cj>° curve, calculated from Eq.(39). 
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FIG. 10. (a) Scaling functions Wh(IC,T), derived from the 
data of Fig5a (p from p = 1CT 1 to 10~ 9 , r < 0.2) by replotting 
in scaling coordinates p 1,/4 /i s //i e vsT = r/p. (b) Kinetic AC— T 
diagram of the regions of power- law asymptotics of Wh(IC, T), 
derived from the data of Fig. 5b (r < 0.2) by replotting in 
scaling coordinates AC = p 3 ^ 2 J vs T = r/p. Note the different 
scale on the T axis in (a) and (b). 



FIG. 11. (a) Scaling functions Wl(IC,T), derived from the 
data of Fig5a (p from p = 10 _1 to 10~ 9 , r < 0.2) by replot- 
ting in scaling coordinates p~ 1 ^ 4 l 3 /l e vs T = r/p. (b) Kinetic 
AC — T diagram of the regions of power-law asymptotics of 
W L {K.,T), derived from the data of Fig.5c (r < 0.2) by re- 
plotting in scaling coordinates AC = p 3 ^ 2 J vs T = r/p. Note 
the different scale on the T axis in (a) and (b) . 



FIG. 12. (a) Scaling functions Wn(IC,T), derived from the 
data of Fig5a (p from p = 10 _1 to 10" 9 , r < 0.2) by re- 
plotting in scaling coordinates p^ 1 ^ 4 N/N e vs T = r/p. (b) 
Kinetic AC— T diagram of the regions of power-law asymptotics 
of W N (K, T), derived from the data of Fig.5c (r < 0.2) by re- 
plotting in scaling coordinates K, = p 3 ^ 2 J vs T = r/p. As 
contrasted from Fig.5c, the regions with n = 1/2 and n = 1 
are added. Note the different scale on the T axis in (a) and 
(b). 



FIG. 13. (a) Collapse of the dependences of hf /h e on J 
from Fig. 6a to the scaling function M(/C) in scaling coor- 
dinates p 1 / A h 1 f /h e vs AC = p 3 / 2 J. Dashed lines show the 
asymptotics (67) and (68). ( b) Collapse of the dependences 
of tm on J from Fig. 6b to the scaling function Tm(AC) in scal- 
ing coordinates tm/p vs AC = p 3 ^ 2 J. Dashed lines show the 
asymptotics (67) and (68). 



FIG. 14. Collapse of calculated numerically for p — 10~ 
(circles) 

and p = 10~ s (diamonds) dependences |$"(J)|,Ti m (J), and 
|mi(J)|, and calculated numerically for J = 10 6 (stars) and 
J = 10 7 (squares) dependences |<^™(p)|, T™(p), and |mi(p)| 
to scaling functions \<f>r(K)\ (a), 7T(AC) (b), and \rrn(K)\ (c), 
respectively Filled symbols stand for i — H, open symbols 
stand for i = L. In Fig. 14a, b filled circle m marks off the 
point AC m = 3.0 with which the AC > AC m region begins where 
a maximum on the \4>l\ curves appears. Open circle c marks 
off the point AC C = 7.3 with which the region AC > AC C of cross- 
ing of the <pL curve with zero begins. Dashed lines show the 
asymptotics (52) and (55) (a), (b), and (54) and (51) (c). 

FIG. 15. Time dependences |<7h(t)| (filled circles) and 
<tl(t) (open circles), calculated numerically at p = 10~ 4 and 
p = 0.5 for J = 10 6 and at p = 0.9 for J = 10 8 (for the 
sake of clarity, the data are given beginning with r = 10~ 2 . 
Straight bold lines show |<7h(t)| and 0l{t), calculated from 
the equations tJi = J 4 i e _ " T , (9), (76), and (74). 



15 



1.6 



1.2 



0.8 



0.4 







